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In complex networks, the failure of one or very few nodes may cause cascading failures. When 
this dynamical process stops in steady state, the size of the giant component formed by remaining 
un-failed nodes can be used to measure the severity of cascading failures, which is critically important 
for estimating the robustness of networks. In this paper, we provide a cascade of overload failure 
model with local load sharing mechanism, and then explore the threshold of node capacity when 
the large-scale cascading failures happen and un-failed nodes in steady state cannot connect to each 
other to form a large connected sub-network. We get the theoretical derivation of this threshold 
in degree-degree uncorrelated networks, and validate the effectiveness of this method in simulation. 
This threshold provide us a guidance to improve the network robustness under the premise of limited 
capacity resource when creating a network and assigning load. Therefore, this threshold is useful and 
important to analyze the robustness of networks. 
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1. Introduction 


Cascading failures are a sort of phenomena that a random failure or intentional attack on 
one or a few nodes lead to serve chain reaction in the networks, which can cause collapse of a 
large fraction of nodes in the network. It is widely found in many real-world networks, such 
as power transmission DEI , communication [3], economic [1], and biological [5j networks. A 
real example of cascading failures is the well-known Northeast Blackout in 2003 [6]. In this 
case, the outage of a generator led to a serve chain reaction of power blackout, which affected 
approximately 50 million people in North America and caused financial losses of about $6 
billion. 

In order to understand the essential mechanism of cascading failures, a number of cascading 
models have been proposed, such as betweenness-based model EH El EE DU, sand-pile model 
mm, and fiber-bundle model JT3i 14], etc. In these models, a node fails when its load 
exceeds its capacity. The failure of this node leads to the redistribution of load in network, and 
can cause collapse of a large fraction of the network. Therefore, the cascading failure process 
of these models is highly depended on the relation between load and capacity. 

When the dynamical process of cascading failures terminates in steady state, the network 
breaks into several connected sub-networks formed by un-failed nodes. The size of the largest 
connected sub-network (i.e., the giant component) can be used to measure the severity of 
cascading failures, which is critically important for estimating the robustness of networks p~5) 
HE] . It is intuitive that when the capacity increases (or the load decreases), the size of cascading 
failures in network reduces, which is confirmed by all of those cascading failure models above. 
In fiber bundle and sandpile model, with the capacity under a critical value (or load above a 
critical value), the giant component disappears, or above which the size of the giant component 
dramatically rises up m H3J- This critical value is a very important feature to measure the 
robustness of networks. However, to our knowledge, there is little research on quantitative 
analysis of the relation between node capacity and this threshold for the break out of cascading 
failures at present. 

In this paper, we provide a cascade of overload failure model with local load sharing mecha¬ 
nism, then explore the threshold of node capacity when the large-scale cascading failures happen 
and there does not exist a large connect sub-network formed by un-failed nodes. 


2. Model 


Here we provide a cascade of overload failure model with local load sharing mechanism. In 
this model, the statuses of nodes are divided into two categories: the un-failed and failed. We 
assume that all nodes in networks are un-failed at the beginning. A node fails if its load exceeds 
its capacity (i.e., overload). When a node fails to work, this node is considered transferring 
a fixed positive load A to each of its un-failed neighbors and being separated from the giant 
component m M- When the cascading failures terminate, only those nodes in the giant 
component are supposed to work. It is natural to assume that the capacity C v of a node v is 
proportional to its initial load U [Z1EI] as 

C v = (1 + oi)L v , (1) 

where the constant a is the tolerance parameter. The initial load of each node is randomly 
distributed following a uniform distribution on the interval [L mini L max \. For simplicity, we 
set L min = 0 and L max = 1 [19] [20] • Then for arbitrary node v, the cumulative distribution 
function of L v is 


0 / < 0 


P{L V <l} = ip{l) = l 0 < / < 1 • 


1 l > 1 


( 2 ) 


The numerical process of cascading failures with local load sharing mechanism is summa¬ 
rized as follows: 


1. Initialization. Generate a degree-degree uncorrelated network with N nodes. Assume all 
nodes are un-failed. The load of each node is uniformly distributed in [L min , L max ], where 
we let L min = 0 and L max = 1. The capacities of nodes are determined by Eq (JTJ). 

2. Beginning. Choose very few node randomly in the network, and set them as failed. 

3. Load redistribution. In each round, each node which is failed in the last round transfers 
a fixed positive load A to each of its un-failed neighbors. An un-failed node turns to be 
failed if it overloads in this round. 




4. Halt. Repeat step [3] if there exist overloaded nodes in the network; otherwise, the process 
halts. Finally, only the un-failed nodes in the giant component are supposed to work. 


3. Analysis 

3.1. Description of critical conditions 

When the cascading failure process ends, the un-failed nodes in the network form several 
connected sub-networks. The fraction of giant component, which is the relative size of the 
largest connected sub-network, can be used as a measure of the network performance against 
cascading failures. Fig. 1 shows the fraction G of the giant component as a function of tolerance 
parameter a in ER random networks [22] and BA scale-free networks |23j . For each G with 
different average degree (k), there exists a critical tolerance parameter a c , with a under which 
the giant component disappears, and over which G increases dramatically and approaches to 
1 finally. When a < a c , G approximately equals to 0, which indicates large-scale cascading 
failures occur and network breaks into extremely small clusters. There are two conditions to 
ensure the G approximately equals to 0 when a < a c : 

I. Large-scale cascading failures occur in the network, i.e., the failed nodes connect to each 
other to form a large connected sub-network. 

II. When the dynamical process of cascading failures terminates in steady state, there does 
not exist a large connected sub-network formed by un-failed nodes. 

In the rest of this section, we give a theoretical derivation of a c with these two conditions 
above. 


3.2. Large-scale cascading failures occur in the network 

First, let us consider condition I that large-scale cascading failures occur in the network. 
This condition equals to the situation that the failed nodes connect to each other to form a large 
connected sub-network. On average, each node fails and causes more than one of its neighbors 
to fail when large-scale cascading failures happen. Consider a node v of degree k v , and any of 
its neighbor w. The probability that node v fails and causes node w to fail is 




(a) ER random networks. (b) BA scale-free networks. 

Fig. 1 The fraction G of giant component changes with increase of tolerance parameter a 

in ER random and BA scale-free networks. 
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where L w is the initial load of node w. 

The cascading failures of nodes are a sort of site percolation process. Because of the locally 
tree-like approximation in percolation of degree-degree uncorrelated networks, the probability 
that node v fails and causes m of its k v — 1 neighbors (subtract the node which causes node v 
to fail) to fail is 


P {Lw + A > (1 + a)L w { — P 



Thus, on average, the node v fails and causes its 



={k v - l)<p 


( 5 ) 


neighbors to fail. 











Denote pk as the probability that a randomly-picked node has degree k. According to Eq ()5]), 
the failure of arbitrary node in the network causes its Y^T=oP k — 1) P (“■) = ((h) — 1) p (A) 
neighbors to fail. When large-scale cascading failures occur in the network, the failure of 
arbitrary node causes more than one of its neighbors to fail, that is 

((k) - l)(p (^j > 1 . (6) 

Thus, for condition 1 that large-scale cascading failures occur in the network, we have 

a <«&>-!) A. (7) 

It is worth noting that there is not any specific assumption on degree distribution for the 
result above. 

3.3. There does not exist a large connected sub-network formed by 
un-failed nodes 

According to the local load sharing mechanism in our model, there may exist some nodes 
who do not propagate cascading failures. These nodes never fail with large tolerance parameter 
a, or their failure will not affect other nodes as all of their neighbors are already failed. We call 
these nodes as absorbing nodes. Before handling with condition II, we pay close attention 
to these absorbing nodes. Let us consider the situation that a node v of degree k v does not fail 
after k v — 1 of its neighbors fail and share loads ( k v — 1)A to it. Consequently, no matter node 
v fails or not, it can not affect other nodes any more. In this situation, node v have the ability 
of absorbing the loads of all its neighbors by itself, thus we call this node as the independent 
absorbing node. The capacity and initial load of an independent absorbing node v meet the 
following relation as 


L v + {k v — 1) A < (1 + o) L v . (8) 

But this is not the only reason that a node happens to be an absorbing node. There is 
another situation for a node to be an absorbing node whose load does not satisfy the Eq (|8j). 
In this situation, considering a node v, m of its neighbors are absorbing nodes, thus we only 
need n v > k v — m — 1 to make node v an absorbing node. We call this kind of nodes as 




(a) The independent absorbing nodes. (b) The dependent absorbing nodes. 

Fig. 2 Two examples of the absorbing nodes. In the procedure of cascading failures, an 

un-failed node (o) may turn to be a failed node (•) when it overloads, or does not fail as it 
is an absorbing node (©). 


dependent absorbing nodes because their abilities of absorbing loads depend on the other 
absorbing nodes of their neighbors. Dependent absorbing nodes can also prevent cascading 
failures locally. 

Fig. 2 gives two examples of absorbing nodes. In Fig. 2(a), the un-failed node 1 has five 
neighbors and four of them (node 2, 3, 4, 5) are failed. Assume that if node 6 does not fail in 
the following procedure of cascading failure, then node 1 does not fail either. Or in another 
case that node 6 fails and transfers a fixed value of load A to node 1. In the latter case, no 
matter node 1 fails or not, it can not affect any other node. Thus, node 1 is an independent 
absorbing node. In Fig. 2(b), we assume that the node 7 is not an independent absorbing node. 
It happens to fail if more than two of its four neighbors fail and transfer loads to it. Assume 
node 8 is an absorbing node, thus it never fail. Consequently, no matter node 7 fails or not, 
it can not share load to any un-failed neighbor. Therefore, node 7 in Fig. 2(b) comes to be 
a dependent absorbing node because its ability of absorbing loads depends on the absorbing 
node 8. 

Then we consider condition II that there does not exist a large connected sub-network 
formed by un-failed nodes. We have explained that absorbing nodes survive and keep un-failed 
after cascading failures. With increase of a, each node has higher probability to be assigned 
with more capacity and be an absorbing node on average. If there are a large number of 
absorbing nodes exist, they may connect to each other to form a large connected sub-network 
which makes G greater than 0. We have divided absorbing nodes into independent absorbing 





nodes and dependent absorbing nodes. Further more, we call an absorbing node as an m- 
absorbing node (m = 0,1, 2, k), if and only if m of its neighbors are absorbing nodes. We 
assume that an independent absorbing nodes is not equivalent to 0-absorbing node. For an 
interdependent node, if it has m absorbing nodes of its neighbors, we still treat this node as an 
m-absorbing node. 

Then, we let the probability that a node of degree k happens to be an absorbing node as 
a*,. The probability that node v happens to be an m-absorbing node is a™' 1 . Then we have 

a k = Yl a fc"° ' ( 9 ) 

m =0 

In a degree-degree uncorrelated network, consider a node v of degree k and any of its 
neighbor w. The probability of node w with degree j is = jp-. Node w happens to be an 

l^i l Pi W 

absorbing node with the probability a y Then, the probability of node n’s arbitrary neighbor 
being an absorbing node is 


O n = 


E 


JPj a j 

(k) 


( 10 ) 


According to Eq. m, the probability that exactly m of node v’s neighbors are absorbing 
nodes is 


( k )<a - ■ ( 11 ) 

\m j 

Node v is an m-absorbing node if it stays un-failed with k — m— loik — m rest neighbors 
fail and transfer loads to it. Therefore, the probability that node v comes to be an m-absorbing 
node with probability 


(m) 

ai = 


<C( 1 - °a) k - m P {L v + (k~m- 1) A < (1 + a)L v } 
k—m I i .. I (k-m- 1)A S 
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( 12 ) 


We substitute Eq. (1T21) into Eq. ((9]), and get the probability that node v (i.e., arbitrary 
node of degree k ) happens to be an absorbing node as 






a«) k - m (l-p{ 


(k — m — 1)A 


a 




elsewhere 


elsewhere 


if k < [ct/A] + 1 


if k < [a/A] + 1 


( 13 ) 


According to the cascade condition [24] [25], the absorbing nodes can not connect to each 


other to form a large connected sub-network with the following condition: 



(14) 


Now we complete the derivation by Eqs. ©, m, ra, and (Thill . These four equations 
depend on each other. We can get the threshold a c with the following iteration process: 

1. Assign (( k) — 1) A as an initial value to a according to Eq. (J7J). 

2. Get a*, via Eq. (fTOl) and (TT3T) with the current value of a. 

3. Substitute ak into Eq. (fT41h If the inequality of Eq. (fT4ll is invalid, discount a very 
small value from a, and then repeat step [21 otherwise, let a c = a and finish the iteration 
process. 

4. Simulation 

We validate the analytical prediction of threshold a c in ER random and BA scale-free 
networks with N = 5000 nodes and varied average degree ( k) = 8,10,12. The load of each 
node is uniformly distributed on the interval [0,1]. We randomly set very few nodes to fail 
as the beginning of the dynamical procedure. When a node fails, it transfers a fixed load 
A = 0.01 to each of its un-failed neighbors. Fig. 3 and Fig. 4 present the giant component 
G with different a, and the threshold a c for G rises from zero to nonzero. Note that different 
values of A only changes the scale of abscissa. For each sub-figure in Fig. 3 and Fig. 4, each 









(a) (fc) = 8 (b) (k) = 10 (c) (k) = 12 

Fig. 3 Validation of the analytical predictions of threshold a c where the giant component 

G exists in ER random networks. With increase of tolerance parameter a , we draw each 

simulated G of 50 realizations respectively. The threshold a c is marked by solid triangle (a). 



(a) (k) = 8 (b) (fc) = 10 (c) (k) = 12 

Fig. 4 Validation of the analytical predictions of threshold a c where the giant component 

G exists in BA scale-free networks. With increase of tolerance parameter a, we draw each 

simulated G of 50 realizations respectively. The threshold a c is marked by solid triangle (a). 

of 50 independent simulation results is draw respectively. The threshold a c is marked by the 
solid triangle. It can be seen that our theoretical predictions are in good agreement with the 
simulation results. 


5. Conclusion 

In this paper, we provide a cascade of overload failure model with local load sharing mech¬ 
anism, and then explore the threshold a c of tolerance parameter for capacity, when the giant 
component G formed by un-failed nodes rises from zero to non-zero. We provide two conditions 
to ensure the G approximately equals to zero when a < a c , which are: I. The large-scale cas¬ 
cading failures occur. II. un-failed nodes in steady state cannot connect to each other to form 






























a large connected sub-network. With these two conditions, we get the theoretical derivation 
of a c in degree-degree uncorrelated networks, and validate the effectiveness of this theoretical 
derivation in simulations. We believe that when creating a network and assigning load with 
local load sharing mechanism, this threshold a c provides us a guidance to improve the network 
robustness under the premise of limited capacity resource. 
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